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C^^ ' We investigate the density structure of the interstellar medium using new 

m ■ 

\0 ' a molecular cloud in the Gum Nebula, made as part of the Spitzer Space Tele- 



O 
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scope Galactic First Look Survey. The maps are correlated with 100 /xm images 

O ■ measured with IRAS. At 24 and 70 /xm, the spatial power spectrum of surface 

brightness follows a power law with spectral index -3.5. At 24 /xm, the power 

law behavior is remarkably consistent from the ~ 0.°2 size of our maps down to 



^ ■ the ~ 5" spatial resolution. Thus, the structure of the 24yum emission is self- 

^ ■ similar even at milliparsec scales. The combined power spectrum produced from 

C^ ■ 

Spitzer 24 /im and IRAS 25 /im images is consistent with a change in the power 

law exponent from -2.6 at spatial wavenumber fc ~ 2 x 10^^ (arcsec)"^ to -3.5 at 
fc ~ 4 X 10~^ (arcsec)"^. The decrease may be due to the transition from a two- 
dimensional to three-dimensional structure. Under this hypothesis, we estimate 
the thickness of the emitting medium to be 0.3 pc. 
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1. Introduction 

More than fifty years ago, astronomers began to reconsider tlie concept of a liomogeneous 
density and velocity structure for the interstellar medium (ISM). The seminal lecture on 
astrophysical turbulence by Chandrasekhar (1949) marked the recognition that studies of the 
ISM were incomplete without knowledge of the physical and observational consequences of 
turbulence. The structure and the dynamics of the ISM affect critically its chemistry (Spaans 
1996; Joulain et al. 1998; Rollig, Hegmann, & Kegel 2002) and its star formation capacity (see 
for example the review by MacLow & Klessen 2004). The ramifications of a turbulent velocity 
and density field for radiative transfer must be taken into account when interpreting spectral 
line and continuum observations of interstellar clouds (Padoan et al. 1998; Hegmann & Kegel 
2000, 2003; Ossenkopf 2002; Juvela & Padoan 2003). Finally, knowledge of the structure of 
Milky Way interstellar matter is essential to the proper interpretation of extragalactic counts 
and measurements of the spatial distribution of the Cosmic Microwave Background (Gautier 
et al. 1992). 

The structural statistics of the ISM are self-similar on a wide range of scales, from 
hundreds of parsecs down to ~ 0.02 pc (Bensch, Stutzki, & Ossenkopf 2001). The angular 
power spectrum of two-dimensional (2-D) images of interstellar emission and absorption 
derived using a variety of tracers yields a power law as a function of wavenumber, k^, with 
exponent ranging from f3 ~ —3.6 (Miville-Deschenes et al. 2003a) to /? ~ —2.5 (Bensch et 
al. 2001) (see also the review by Falgarone, Hily-Blant, & Levrier 2004). Power spectral 
analysis of the far-infrared emission from the Galactic cirrus results in the same range of 
values on angular scales of 1' and larger (Gautier et al. 1992; Herbstmeier et al. 1998). 

The Spitzer Space Telescope allows us to extend the examination of turbulent density 
fields down to ~ 5" scales. We present new high-resolution maps of the 8 /im, 24 /im, and 
70 /xm diffuse emission towards the Gum Nebula, made as part of the Spitzer Galactic First 
Look Survey (§2). In §3 we derive colors between the Spitzer and IRAS bands and perform a 
power spectral analysis of the Spitzer images. For the 24 /im and 70 /im maps we find power 
law exponents similar to those derived from H I observations, with no detectable break at 
high wavenumber down to spatial scales of ~ 0.01 pc. We discuss implications of our results 
for interstellar structure in 54. 
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2. Observations 

2.1. The Spitzer Galactic First Look Survey 

The new data that we analyze form part of the Spitzer Space Telescope (Werner et al. 
2004) Galactic First Look Survey (GFLS) (Noriega-Crespo et al. 2004)^. We describe here ob- 
servations made using the IRAC (Fazio et al. 2004) and MIPS (Rieke et al. 2004) instruments 
on board Spitzer. The IRAC observations (PID 104; ads/sa.spitzer^6579712) were con- 
ducted on 2003 December 7, and the MIPS observations (PID 104; ads/sa.spitzer7^6578176) 
were conducted on 2003 December 9. 

The surface brightness maps on which we have performed power spectral analysis are 
shown in Figure 1. The images on the left side of the figure are the output of the automated 
post-Basic Calibrated Data (post-BCD) Spitzer calibration pipeline^. For each of the three 
Spitzer bands considered in this study — IRAC 8/im and MIPS 24 and 70 /xm — we have ana- 
lyzed slightly different fields, albeit with considerable overlap. There was a small mismatch 
in the observed IRAC and MIPS fields caused by a ~ 2° rotation of Spitzer between the two 
observation dates. We also truncated the 70 /xm mosaic due to excessive noise in portions of 
the image. 

Each of the examined regions are centered on Galactic coordinates (/, b) = (254.5, —9.5) 
and cover a square of size ~ 0.°1 — 0.°3 (see Figure 1), comprising less than 2% of the GFLS 
total sky coverage. This line of sight intersects a molecular cloud associated with the Gum 
Nebula, an expanding supershell of radius ~ 70 — 130 pc (Yamaguchi et al. 1999). The CO 
gas associated with the Spitzer field has negative radial velocities, placing it on the near 
side of the bubble. Since the expansion center is at most 500 pc from the Sun (Woermann, 
Gaylard, & Otrupcek 2001), we estimate a distance d ~ 400 pc to our field. The 24/i,m and 
70 /im fields also include the object DC 254.5-9.6 from the optical catalog of southern dark 
clouds of Hartley et al. (1986). 

Before estimating the power spectrum, the Spitzer images shown in the left column 
of Figure 1 were processed to remove point sources and instrumental noise. The emission 
from point sources was characterized and removed using the StarFinder code written in IDL 
(Diolaiti et al. 2000a,b). The images were then filtered to remove instrumental noise using 
a multiresolution wavelet technique (Starck & Murtagh 1998). The method required an 
estimate of the noise spatial response. Following Miville-Deschenes, Levrier, & Falgarone 



^http://ssc. spitzer. caltech.edu/fls/galac/ 
^http://ssc. spitzer. caltech.edu/postbcd 
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(2003b), we used the uncertainty images automatically produced by the post-BCD pipeline 
to account for spatial variations in the noise. The right column of Figure 1 shows the left 
column images after point source removal and noise filtering. 

Power spectra for each of the images in Figure 1 were computed by squaring the am- 
plitudes of the image fourier transforms. Prior to transforming the images, the outer 5% 
of pixels were multiplied by a cosine taper. This apodization function minimized edge ef- 
fects caused by our non-periodic data. In fourier wavenumber (k) space, the effect is a 
slight smoothing of the amplitudes. The radial power spectra were derived by azimuthally 
averaging the squared fourier amplitudes in equally-spaced wavenumber bins. Statistical un- 
certainties were calculated as Xj ^fnl times the azimuthal average, where n^ is the number of 
amplitude samples comprising the average. To ensure that we were not adding structure by 
noise filtering, the insignificant wavelet coefficients and scales that had been rejected in the 
filtering process were summed to make a "noise image." Power spectra of the noise images 
showed no systematic trends as a function of spatial wavenumber, consistent with uncor- 
related "white" noise. In addition, a subregion of the Spitzer mosaics with no detectable 
emission showed the same power spectrum as the noise images. 



2.2. IRAS Sky Survey Atlas Images 

We used surface brightness maps at 12, 25, 60, and 100 /xm from the IRAS Sky Survey 
Atlas (ISSA). ISSA maps were produced that covered a 12° x 12° region centered on the 
Spitzer subregions. These maps were corrected for gain and offset variations by comparing 
with measurements taken with the Cosmic Background Explorer (COBE) Diffuse Infrared 
Background Experiment (DIRBE) (Miville-Deschenes & Lagache 2004). Before power spec- 
tral estimation, the IRAS data were processed for source removal, noise filtering, and PSF 
deconvolution according to the method described by Miville-Deschenes, Lagache, & Puget 
(2002). 



3. Results 

3.1. Spitzer and IRAS colors 

Table 1 lists mean surface brightness ratios Ix/hoo between the Spitzer and IRAS data. 
The point source-subtracted Spitzer maps were convolved with the IRAS beam and resam- 
pled onto the ISSA grid. Linear fits of I\ as a function of Jioo were performed. The slopes 
and their errors are reported in Table 1. For most wavelengths the data were well correlated 
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with Jioo- After processing, the MIPS 70 /xm image (Figure 1 [f]) still showed systematic 
instrumental features, most notably vertical stripes in the direction of the scan legs. We do 
not analyze the 70/100 color here, in expectation of improved processing of the 70 /xm data. 
The Spitzer-derived spectral energy distribution of the cirrus will be discussed more fully in 
a future paper. 



3.2. Spatial Power Spectrum 

Figure 2 displays the spatial power spectra of emission in the fields shown in Figure 
1. In each of the three Spitzer data sets (IRAC 8/im, MIPS 24 and 70 /xm) point source 
subtraction and noise filtering had an obvious effect on the power spectra. Both processing 
steps removed a mostly flat "white" noise component. The benefits of point source removal 
were most pronounced in the IRAC 8/xm data, while noise filtering improved greatly the 
linearity of the spectra (in log-log plots) in all three fields. 

Power law fits (P = Ak^^) to the point source removed/noise filtered spectra are super- 
imposed on the data in Figure 2. The 24 /im and 70 /im fits gave the same result within the 
formal fit errors, (^24 = -(3.52 ± 0.01) and pro = -(3.48 ± 0.04). In log-log plots, the 24 /xm 
spectrum did not deviate from a straight line for almost two decades of spatial frequency, 
down to the spatial resolution of ~ 5". At the ~ 400 pc distance of the emitting material, 
this is equivalent to a size scale of ~ 0.01 pc. 

In contrast, the 8;um image had a power spectrum over the same spatial scales that was 
flatter by 0.4, [3^ = —(3.13 ±0.03), indicating more power at small scales than the other two 
images. A bump in the 8 /im power spectrum at high wavenumber was probably caused 
by incomplete source subtraction. Some sources were clearly not pointlike and subtracting 
a scaled point spread function added high frequency structure. Unfortunately, this occured 
where the IRAC 8/im image might have provided additional information on the structure 
at the smallest scales. 

We extended the power spectrum of our 24 /im field to large scales using the power 
spectrum of the IRAS map that includes it. Figure 3 shows a comparison of the MIPS 
24 /xm power spectrum with that of an ISSA 25 /xm image that includes the same field. The 
25 /xm data were normalized (multiplied by a constant) to match the 24 /xm spectrum in the 
overlap region and all data points were rebinned for clarity (except for the low-A; data, each 
plot symbol is an average over more than one radial point). The 25;um power spectrum was 
fit by a power law with j325 = —(2.63 ± 0.02), which is shallower than the 24 /xm spectrum 
by ~ 1. Our data suggest that there is a transition in the power law exponent from -2.6 
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to -3.5 somewhere in the wavenumber range 2 x 10"^ < A; < 4 x 10~^ (arcsec)"^. This also 
seems to occur in the 70 /xm data. If the decrease in slope is real, then it is comparable to 
the bend in the power spectrum of H I seen by Elmegreen, Kim, & Stavely-Smith (2001) 
in the Large Magellanic Cloud (LMC), where the spectral exponent also decreased by ~ 1 
from low to high k. This is predicted to occur when the transverse size of an image is 
greater than the line of sight depth of the emitting medium, ^los- More precisely, the 
transition occurs at wavenumber k = l/2di^os (for theoretical and numerical support for 
this argument, see Lazarian & Pogosyan 2000; Miville-Deschenes et al. 2003b). Under this 
hypothesis, we estimated c^los for the 25 and 24 /xm-emitting dust cloud. Assume that the 
transition occurs at /Ctrans = 0.003 (arcsec)"^, or 1.54 pc~^ if the distance is 400 pc. Thus, 
c/los = l/2/i;trans = 0.3 pc. Wlthout invoking distance we can use the angular size of our map, 
©map; ciud the trausitlou wavenumber to estimate the size-to-thickness ratio of the medium, 
L/di^os = 26inap^trans- Taking 6map = 0.°3, i.e., the size of the 24 /xm image used to compute 
the high-/c power spectrum, L/di^os = 6.5. This is probably a lower limit, since the emission 
continues beyond the edges of the map. 



4. Discussion 

The Spitzer 24 /xm and 70 /im surface brightness maps of the Gum Nebula have the 
same power spectral index, (3 ~ —3.5. The 24yum data in particular conform well to a -3.5 
power law down to the spatial resolution of ~ 5", equivalent to a size scale of ~ 0.01 pc at 
the estimated distance of ~ 400 pc. Thus, the self-similar structure of the 24ij-t^ emission 
continues to milliparsec scales. To understand properly the relationship between the emission 
and the column density, we need to model the variations in grain heating and emissivity as 
a function of density, extinction, and distance from ultraviolet sources. Nevertheless, these 
f3 estimates that we have derived for the mid-infrared dust surface brightness are close to 
the value —3.6 reported for the density and velocity fields of H I in emission in the Galaxy 
and the LMC (see Miville-Deschenes et al. 2003a; Elmegreen et al. 2001, respectively). A 
spatial power spectrum index of —3.5 is also close to the prediction P = —11/3, based on the 
incompressible turbulence theory of Kolmogorov (1941). In contrast, the velocity integrated 
emission from dense molecular gas typically yields values of /9 ~ —2.8 (Stutzki et al. 1998; 
Bensch et al. 2001; Padoan et al. 2004). 

Our power spectral analysis may have enabled us to access the depth of 24 /xm emitting 
material along the line of sight, c/los ~ 0.3 pc. The steepening of the power spectrum from 
/92D to PsB ~ P2B — 1 (in our case from -2.6 to -3.5) is predicted to occur when the size of the 
map is larger than di^os (Miville-Deschenes et al. 2003b). At small wavenumbers the image 
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power spectrum reverts to that of a 2-D field, i.e., the projected 3-D volume statistics are 
equivalent to those of a 2-D slice, and the power spectrum flattens.^ The 2D/3D transition 
has been observed before in the LMC (Elmegreen et al. 2001), but this is the first example 
in the Milky Way. 

What is the physical significance of the 0.3 pc emitting layer thickness? There are two 
possibilities that depend on the excitation of dust grains into emission in the 24 /im band. If 
the emission is dominated by dark regions, then c^los could be the average "skin depth" for 
penetration of UV photons, or the depth beyond which the abundance of emitters decreases 
drastically. The skin depth effect should be even more apparent in the 8/xm emission, but 
less apparent in the 70 /xm emission. 

If, on the other hand, the emission comes from mostly diffuse regions with A^ < 1 mag, 
then UV photons permeate the entire cloud and c^los is the average thickness of the Gum 
Nebula shell itself. While opaque lines of sight do exist in the Gum Nebula (e.g., the 
cometary globules studied by Sridharan 1992), on ~ 2' scales our field has only moderate 
extinction, A^ < 1.1 mag (derived from the brightest ^^CO contour of Yamaguchi et al. 
1999). Furthermore, the 70;um power spectrum also shows evidence of a break at /c = 
0.003 (arcsec)^^, implying that c^los = 0.3 pc is independent of excitation. The 24 and 70um 
data show only one scale where the 2D-3D transition occurs. 

It would not be surprising if the 24 /im size-to-thickness ratio, L/c^los ^ 6.5, applied 
not only to the small portion of the Gum Nebula that we have observed, but to the diffuse 
medium of the Galaxy in general. In an H I survey of the cold neutral medium, Heiles & 
Troland (2003) demonstrated that clouds are "sheetlike," with size-to-thickness ratios of up 
to 280. Far from being unusual, the 2D/3D power spectrum break that we see may be a 
normal feature of large spatial dynamic range observations of the local cold neutral ISM. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. We thank the referee, Paolo Padoan, for his comments that improved the 
manuscript. 



•^An alternative interpretation holds that the break in the spectral index is caused by the lack of self- 
similarity in the 3-D turbulence above the scale where the system becomes two-dimensional (Padoan et al. 
2001). Even under this interpretation, however, the observed break would mark the thickness of the system. 
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Fig. 1. — Spitzer surface brightness maps used for power spectral analysis. The field includes 
a molecular cloud associated with the Gum Nebula, {(a) and (b)): IRAC 8/xm gray levels 
range from 0.28 (white) to 0.33 MJysr"^ (black), {(c) and (d)): MIPS 24 /xm gray levels 
range from 19.1 to 20.8 MJysr"^. {(e) and (f)): MIPS 70 /xm gray levels range from 20 to 
490 MJysr^^. The left column — frames (a), (c), and (e) — depicts the original post-BCD 
mosaic images. The right column — frames (b), (d), and (f) — shows the results of point source 
subtraction and noise filtering (see text) applied to the left column images. Each map has 
the same grid of Galactic coordinates superimposed (grid lines are separated by 0.1 degree). 
A cross in each of (c) through (f) shows the position of cloud DC 254.5-9.6 (Hartley et al. 
1986). 
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Fig. 2. — Spatial power spectra of emission towards the Spitzer fields in Figure 1: {top) IRAC 
8/im; {middle) MIPS 24 /im; {bottom) MIPS 70 /im. Spectra correspond to the two stages of 
image processing depicted in Figure 1: original post-BCD mosaic (dotted lines), and point 
source removed/noise-filtered (solid lines). Power law fits to the point source removed/noise- 
filtered data are superimposed as gray lines and labeled on the graphs. Statistical error bars 
are plotted for the point source removed/noise-filtered measurements. 
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Fig. 3. — Spatial power spectra at 25;um {open triangles; IRAS) and 24yum {filled squares; 
Spitzer) for the field encompassing the images displayed in Figure 1. Both spectra have been 
point source- removed and noise- filtered, albeit with different methods (see §2). The 25 /im 
spectrum has been normalized to intersect the 24 fira spectrum. The data have been rebinned 
in equally-spaced logarithmic intervals for clarity. Separate power law fits are superimposed 
as dashed lines on the two spectra, and labeled on the graph. We do not plot statistical 
error bars since they would be smaller than the symbols. 
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Table 1: Colors in Spitzer and IRAS Maps 



A 


h/hoo^ 


(/im) 




8 


0.039 ±0.005 


12 


0.029 ±0.005 


24 


0.041 ±0.005 


25 


0.048 ± 0.003 


60 


0.315 ±0.003 



"Slope of linear fit to the surface brightness I\ as a function of /loo, ± 
formal fit errors. 



